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The magnetotransport in single layer graphene has been experimentally investigated in magnetic 
fields up to 18 T as a function of temperature. A pronounced T-dependence is observed for T < 50 
K, which is either metallic, or insulating, depending on the filling factor v. The metal-insulator 
transition (MIT) occurs at \v c \ ~ 0.65 and in the regime of the dissipative transport, where the 
longitudinal resistance R xx > \Rk- The critical resistivity (R xx per square) is p xx (v c ) ~ \Rk and 
is correlated with the appearance of zero plateau in Hall conductivity <j xy (y) and peaks in a xx (v). 
This leads us to construct a universal low-T (n, B) phase diagram of this quantum phase transition. 
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Graphene, a honeycomb crystalline single layer of car- 
bon atoms, is a gapless two-dimensional (2D) semimetal 
with many novel electronic properties [lH^]. The low- 
energy charge excitations in graphene can be described 
by the Dirac equation for massless chiral fermions. In 
a perpendicular magnetic field B the electronic energy 
eigenstates become quantized Landau levels (LL) whose 
energies are given by En = sgn(N)hw c y/\N\, where 
w c = VF\/2eB /(he) is a "cyclotron frequency", vf is 
the fermion velocity and N = 0, ±1, ±2, ±3... is an inte- 
ger. Unlike in a non-relativistic 2D electron gas (2DEG) 
system, there is a field-independent LL at E = for 
N = in graphene, which has a four-fold degeneracy, 
corresponding to two Fermi points of the honeycomb lat- 
tice (K and K 1 valleys), and two spin states. The N = 
level is a fundamental property of the relativistic Dirac 
Hamiltonian and is particle-hole symmetric. As a result, 
an anomalous quantum Hall effect (QHE) has been ob- 
served in graphene, whose quantized Hall resistance is 
given by R xy = \RkI{N +5), where Rk = h/e 2 is the 
resistance quantum [2J, |3j. A plateau at R xy = Rk/2 ap- 
pears at a filling factor v = n s h/(eB) — ±2, where n s is 
the charge density per unit area, because only two of the 
four N = LL are available to either electrons, or holes. 
In the semiclassical description of magnetotransport this 
leads to a Berry phase 7r, characteristic of Dirac particles. 

Electronic states with different spin and valley indices 
split even further (see Fig. Q](a)). Spin-splitting results 
from the Zeeman interaction, and K and K' valley split- 
ting from lattice imperfections and inter-valley scatter- 
ing inherent to the mapping of the lattice Hamiltonian 
of graphene onto axially symmetric Dirac fermions 0. 
The nature of the electronic states at the N = LL has 
recently been under intense theoretical and experimental 
scrutiny, but it still remains unclear. While some theo- 
ries predict unusual metallic transport via gapless edge 
states (fUll other suggest that a gap may open at high 
magnetic fields, resulting in an insulating state near the 
charge-neutrality point (CNP), where n s — 0- 11 1. 



Several groups have studied experimentally charge 
transport in the N — LL regime near the CNP in 
graphene, with conflicting results [21 l3l la. Il2l4l7l |. Some 
groups have observed finite longitudinal resistance R xx < 
Rk even at high magnetic fields, where the QHE plateau 
at R xy = Rk/2 is well developed [2J, [j, |8|, [l3|, ]l4| . In con- 
trast, others have found a marked resistance increase, to 
R xx 3> Rk, and a plateau of nearly zero Hall conductiv- 
ity near the CNP in high fields, suggesting a breakdown 
of the N = QHE and an insulating state at low tem- 
perature [Til. [T5HT7I . 

In some studies a plateau at v = 1 and R xy = Rk was 
also observed, indicating inter-valley splitting of the N = 
LL, [12]. In Ref. [15], R XX {B) > 100R K was reported 
at the CNP, and its divergence with field was associated 
with a Kosterlitz-Thouless (KT) transition at a sample- 
dependent critical field B c (> 18 T), which was lower 
for better quality samples. The field-induced insulator 
behavior was later correlated with the filling factor v, 
and evidence for two different insula ting phases emerging 
with decreasing v < 1 was found [17| . More recently, 
magnetic field driven insulator phases have been observed 
in both single and double layer high mobility suspended 
graphene samples, where they emerge from v = 1, as well 
as fractional, v = 1/3 QHE states [18h20I |. 

The phase diagram of the electronic states at the 
N = LL in graphene emerging from the above obser- 
vations is quite complex. However, it seems to have a 
universal constituent, a metal-insulator transition (MIT) 
where the incompressible electronic liquid state of the 
QHE plateau is followed by a localized insulating phase 
[Fig. Q] (b)]. How does dissipative regime develop in 
a system of Dirac electrons in graphene, where robust 
quantum Hall transport is sustained up to hunderd K and 
above, is a fundamental open issue, which we address here 
by performing temperature-dependent magnetotransport 
measurements. We identify the first of these MIT, which 
occurs at v w 0.65, and establish its low-temperature (n s , 
B) phase diagram, shown in Fig. Q] (c). 
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Figure 1: (a) Splitting of the Landau levels in graphene. (b) 
Our proposed schematic phase diagram of the electronic states 
on the N = LL consistent with the results of Refs. [l2l | - 
[20l |. (c) The experimental phase diagram of the Quantum 
Hall Metal-Insulator transition in graphene obtained from our 
data. It also accommodates various data published in Refs. 
03 - HU, reconciling existing observations. Symbols are ex- 
plained in the text. 

We have studied graphene samples prepared by me- 
chanical exfoliation of pyrolytic graphite [lj deposited 
onto a highly doped silicon substrate covered with 285 nm 
thick SiC>2 layer. Hall bar devices with Cr(3 nm)/Au(50 
nm) contacts were patterned and etched using electron 
beam lithography and oxygen gas plasma etching, with 
the Si substrate serving as a gate to control the nature 
(electrons or holes) and density of the electronic carri- 
ers. (The optical image of the device investigated in this 
study is shown Fig. [3] (a)). Prior to measurements, de- 
vices were annealed for several hours in vacuum at 400 
K, to remove impurities. Four-probe magnetoresistance 
measurements were performed using a low frequency (7 
Hz) lock-in technique with a / — 10 nA current. 

The longitudinal and Hall resistivities p xx t xy > |2l| of 
our graphene sample at the base temperature (T = 0.25 
K) and B = 18 T are shown in Fig. [2] The correspond- 
ing conductivities were obtained by standard matrix in- 
version, a xx{yx) = p xx ( xy )/{p 2 xx + ply)- A peculiar be- 
havior is observed near N = LL filling, where p xx (and 
R xx ps 3p xx ) become large, noticeably exceeding Rk, 
and a two-peak structure in a xx concomitant with a "zero 
plateau" of a xy ~ appears [Fig. [2] (b)]. These observa- 
tions agree with previously reported findings [i , 15 - 17 1 . 

The full (n s , B) dependence of the longitudinal resis- 
tivity and conductivity at T = 0.25 K is summarized in 
Fig. [3J in the form of color contour plots. The data was 
collected by gate voltage sweeps at fixed B every 2 T. 
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Figure 2: QHE data as a function of the gate voltage V g , for 
B = 18 T at T = 0.25 K, used in Fig. El (a) p xx and p xy , top 
axis shows the corresponding carrier density n s ; (b) a xx and 
<j xy , top axis shows LL filling factor v. 

To avoid errors from contacts misalignment, p xy and p xx 
were measured for two opposite directions of magnetic 
field and then antisymmetrized and symmetrized, respec- 
tively. As seen in Fig. [3] (a), p xx displays well-developed 
Shubnikov-de Haas (SdH) oscillations fanning away from 
the CNP, Vo ps 7.5 V, with increasing B. Maxima of p xx 
track LL filling v — 0, ±2, ±6, ±10, which is typical 
of graphene Despite the moderate mobility of our 

sample, (/i ps 5.8- 10 3 cm 2 /Vs), the QHE regime develops 
for B > 4 T. with plateaus at p xy ps \Rk and p xx ps 0. 

A striking feature in Fig. [3J (a) is a well-defined resis- 
tive region near n s = that appears at B > 9 T, where 
Pxx ^ Rk/2- Its boundary is clearly correlated with 
maxima in a xx (n s , B) [Fig.[3J(b)], and in the derivative of 
the Hall conductivity, da xy (ns, B)/dn s [Fig. [3J (c)]. The 
latter implies step- like behavior of a xy , corresponding to 
a plateau phase near n s = 0. While these observations 
clearly identify a well-defined region in (n s , B) space cor- 
responding to a resistive plateau phase near n s = for B 
> 9 T, its nature remains unclear. Is this indeed a new 
electronic phase (insulator?) at the N = LL, bounded 
by a phase-transition line in the (n s , B) phase diagram? 

In order to answer this question and understand the 
appearance of the resistive N = LL state in Fig. [3J we 
studied the temperature dependence of the magnetore- 
sistance data at B = 18 T shown in Fig. O The results 
of gate voltage sweeps at T = 0.5, 1.5, 4, 10, 20, 30 and 
50 K are shown in Fig. 2J Two regions are easily iden- 
tified in the figure: a low-n s region where the longitudi- 
nal resistance increases with decreasing temperature, and 
another region at high carrier densities in which the resis- 
tance shows the opposite, metallic behavior [Fig. |U(a)]. 
An obvious feature separating the two regions are cross- 
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ing points of the curves measured at different tempera- 
tures, where p xx is temperature-independent. Inspection 
of Fig. [5] (b) , (c) shows that the insulator phase at low 
n s is marked by the appearance of a a xy ~ plateau 
and an accompanying two-peak structure in a xx . Thus, 
a low-ro s plateau phase identified by the behavior of the 
conductivity in Figures [2 [3] at T = 0.25 K for B > 9 T, 
indeed corresponds to a distinct, insulating phase. 

Metal-insulator transition in Fig. 0] occurs for vf^ s» 



0.57, 



0.5Rk and v, 
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showing some 20% asymmetry between the electron and 
hole sides, respectively. However, the critical resistivity 
is remarkably close to the value Rk/2 expected for a 
metal-insulator transition in disordered 2D electron sys- 
tems (the denominator 2 is a consequence of graphene 
valley degeneracy) [22j . Since contacts in our etched Hall 
bar device are non-invasive, the observed electron-hole 
asymmetry is unlikely to be contact-related. Perhaps, 
it arises from breaking of the particle-hole symmetry by 
random external potential, suggesting a connection of the 
observed MIT with electrostatic disorder. 

A similar behavior of resistivity isotherms measured 
at different T function of magnetic field has pre- 
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Figure 3: (a) Contour map of the longitudinal resistivity, 
Pxx = Rxx ■ Ws/L s , as a function of charge density n 3 and 
magnetic field B, L s — f.5(f)/j,m, W 3 = 0.5(f)/im are sam- 
ple dimensions. Top axis shows the back gate voltage V g . 
LL are clearly identified by ridges fanning away from n s = 
with increasing B, traced by black broken lines. Left bottom 
inset is the optical image of the device. Right bottom inset 
shows the LL fan diagram with fits yielding n 3 = Cv(V g — Vo), 
Cv = 9.2(4) ■ f 10 cm _2 /V, for charge density and Vo = 7.5(f) 
V for the charge neutrality point (CNP). (b) The correspond- 
ing maps of the longitudinal conductivity, a xx , and (c) the 
derivative of the Hall conductivity, daxy/dn s . 




Figure 4: Temperature dependence of the longitudinal resis- 
tivity (a), and longitudinal (b), and Hall (c), conductivities 
versus V g and carrier density (top axes) at B = f 8 T. Vertical 
broken lines show critical densities separating the insulating 
and the metallic phases, and are determined from the crossing 
points in panel (a). 

viously been observed in conventional 2DEG systems 
at low carrier densities 23- 27| • The crossing point of 
the p xx {B,T) curves where the sign of the temperature 
derivative dp xx /dT changes, has been associated with a 
transition from the electronic liquid state of the QHE 
plateau phase to an insulator, as the LL filling is de- 
pleted with increasing B. Such plateau-insulator tran- 
sitions (PIT) have been observed both for integer and 
fractional QHE plateaus and are understood as metal- 
insulator quantum phase transition induced by Ander- 
son (strong) localization of electrons. In conventional 
2DEGs the transition from the last integer QHE plateau 
to an insulator occurs at p xx ~ Rk and in the range 
0.52 < v c < 0.8, similar to what we find in graphene, 
even though massless charge carriers are expected to 
avoid strong localization and the localization mechanism 
in graphene must be of different nature. 

Association of the MIT with the crossing point of resis- 
tivity isotherms is only meaningful when p xx (B, T) shows 
clear metallic and insulating behavior in a finite temper- 
ature range on the two sides of the transition (2?]]. As 
shown in Fig. [5J this is indeed the case for 1 K < T < 50 
K for the B = 18 T data in our graphene sample. In a 
temperature window 10 K < T < 50 K, log p xx (B, T) is 
consistent with nearly linear dependence on T -1 / 2 [Fig. 
[5] (b)], indicative of Coulomb-gap-induced variable range 
hopping (VRH). Despite the limited dynamical range of 
our data, such dependence agrees visibly better than a 
simple activation dependence T" 1 . This provides a hint 
that scattering mechanism driving the insulator state is 
a screened Coulomb interaction. Remarkably, the same 
T-dependence is retained across the transition, persisting 
with the slope sign change on the metallic side. Similar 
symmetry, where p xx {T,v - v c ) oc l/p xx (T v r - v), has 
also been observed in low mobility 2DEG [25| . 

Having established the signature of the MIT on the 
N = LL in graphene, we now return to construct- 
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Figure 5: (a) Semi-log plot of p xx versus temperature at 
B — 18 T for 12 filling factors i/, revealing the metal-insulator 
transition at v c ~ —0.7 and 0.6, where p xx {v c ) ~ Rk /2 
(dashed line), (b) Same data as a function of T~ 1//2 . 

ing its low-T (B, n s ) phase diagram. This is shown in 
Fig. Q] (c), which complies various data reported ear- 
lier by other groups along our own results. Closed sym- 
bols show transition identified from the longitudinal re- 
sistivity, using p xx — Rk/2, while open symbols cor- 
respond to the <T xy ~ plateau phase, identified from 
o'xyi^c) — i- e. the mid-point between the 

N = and N — 1 plateaus. The overall topology of the 
phase diagram, including the electron-hole asymmetry, 
is remarkably universal, but the exact MIT critical fields 
are slightly sample-dependent, especially at low B. 

In conclusion, we have presented evidence for a metal- 
insulator transition on the N = LL in graphene, which 
occurs in the regime of the dissipative transport, where 
Rxx > Rk/2, and at p xx « R K /2 (R xx ks 3p xx for our 
sample). It is surprisi ngly similar to the plateau-insulator 
transition in 2DEG 231427 1. even though strong local- 
ization driving the PIT in 2DEGs should be absent in 
graphene. This difference is clearly seen in the low-field 
part of the phase diagram of Fig. [TJ where Klein tun- 
neling of Dirac electrons in graphene impedes the low- 
density insulating phases [27j. In 2DEGs, the conductiv- 
ities in the vicinity of the PIT follow a semicircle relation 
[28| . indicating a transition to a peculiar quantum Hall 
insulator state [25|, [2j| . It was found that in graphene 
a semicircle describes very well the transitions between 
plateaus in the metallic phase For the N = LL it 

predicts p xx +p xy — (Rk/2) 2 , and therefore p xx < Rk/2- 
We find that this is violated in the insulating phase, sug- 
gesting that p xx m Rk /2 (a xx ~ 2e 2 /h) is the maximum 
metallic resistivity (minimum conductivity) of the dissi- 
pative Hall state on the N — LL in graphene. 
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